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FOREWORD 


This  technical  report  was  prepared  for  the  Director 
of  Aerospace  Sciences,  Air  Force  Office  of  Scientific 
Research,  Washington,  D.  C. ,  in  completion  of  Con¬ 
tract  AF49(638)-1575.  It  describes  research  studies 
carried  cut  by  the  fttysical  Processe  Group  of  the 
Research  Division  of  Rocketdyne,  a  Division  of  North 
American  Aviation,  Tne. 
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ABSTBACT 

An  analytical  and  experimental  program  to  investigate  the  effects  of 
aluminum  on  solid-propellant  combustion  instability  was  conducted.  The 
analytical  phase  was  devoted  to  a  determination  of  the  energy  dissipation 
rate  associated  with  possible  damping  mechanisms  at  either  the  burning 
propellant  surface  or  in  the  chamber  volume.  Particulate  damping,  droplet 
deformation,  surface  modes  of  vibration,  the  effects  of  aluminum  on  the 
shattering;  and  combustion  interaction  of  the  burning  surface  with  the 
acoustic  field  were  all  considered.  The  relative  magnitude  of  these 
damping  mechanisms  are  developed  and  discussed. 

The  effects  of  aluminum  on  nonacoustic  instability  were  also  considered 
in  the  analytical  phase.  It  was  shown  that  the  nonacoustic  problem  could 
be  treated  by  employing  conventional  acoustic  methods,  token  to  an  appro¬ 
priate  limit.  This  treatment  was  applied  to  various  existing  experimental 
data  on  nonacoustic  instability,  and  good  agreement  was  obtained.  This 
information  was  used  to  propise  a  mechanism  for  quasi -periodic  shedding  of 
aluminum  from  the  burning  surface. 

This  analysis  was  followed  by  an  experimental  phase  having  the  goal  of 
separating  energy  dissipation  mechanisms  into  surface  and  chamber-volume 
effects  and  establishing  the  relative  magnitudes  of  each.  These  measure¬ 
ments  were  made  employing  a  T-burner  in  which  a  pressure  pulse  could  be 
generated.  Measurements  were  made  on  propellants  containing  up  to  15-percent 
aluminum  at  mean-pressure  levels  up  to  600  psi  and  over  a  frequency  range  of 
200  to  2000  cps.  A  discussion  is  presented  of  the  results  of  these  experi¬ 
ments,  particularly  with  respect  to  the  damping  mechanisms  considered  in  the 
analysis.  Some  information  applicable  to  nonacoustic  instability  was 
obtained  as  well. 
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INTRODUCTION  AND  SUMMABY 


Combustion  instability  has  persisted  as  g  problem  in  motor  development 
programs  for  many  years.  The  addition  of  powdered  metals  has  been  found 
to  be  highly  effective  in  suppressing  acoustic,  or  high-frequency,  instabil¬ 
ity.  Aluminum,  in  particular,  is  the  principle  means  of  suppressing  this 
form  of  instability,  but  the  presence  of  aluminum  appears  to  intensify 
nonacoustic-instability  problem.  In  spite  of  its  extensive  use,  the 
mechanisms  by  which  aluminum  operates  as  a  stabilizing  influence  are  not 
well  understood.  (Ref.  1  through  3). 

One  mechanism  suggested  for  the  stabilizing  influence  of  aluminum  is  the 
damping  associated  with  the  condensed  phase  products  of  the  aluminum 
combustion.  The  knowai  effectiveness  of  particulate  damping  indicates 
that  this  mechanism  must  play  an  important,  if  not  dominant ,  role,  also, 
interaction  of  the  aluminum  with  the  combustion  has  been  suggested  as 
a  stabilizing  influence  (Ref.  i  through  8).  Experimental  evidence  has 
been  cited  for  both  of  these  effects  but,  except  in  one  instance,  no 
definitive  answers  have  evolved. 

For  low  aluminum  concentrations,  at  least  for  one  propellant  system  and 
for  up  to  1.5  weight-percent  of  aluminum,  particulate  damping  has  been 
found  to  be  dominant  (Ref . 0  and  7).  However,  practical  propellants 
frequently  contain  10-  to  20-percent  aluminum.  At  these  loadings  the 
balance  between  particulate  damping  and  other  damping  mechanisms  may  be 
completely  shifted.  To  use  powdered  metals  most  effectively  to  promote 
stability,  the  stabilizing  processes  must  be  better  understood.  •  With 
this  motivation,  a  research  program  was  undertaken  to  investigate  the 
effect  of  aluminization  on  solid  propellant  combustion,  and  in  particular, 
to  differentiate  between  surface  and  chamber-volume  effects  in  damping  cr 
driving  combustion  oscillations. 
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The  initial  phases  of  tin*  program  wort*  theoretical.  Various  processt-s 
wore  investigated  in  an  effort  to  estimate  analytically  their  relative 
importance,  witltin  each  of  the  categories  cited  (surface  and  chamber- 
volume  effects),  as  well  as  with  respect  to  each  other.  In  some  cases 
these  estimates  could  be  made  on  the  basis  of  a  fairly  well-defined 
theory.  In  others  the  complexity  and  lack  of  understanding  of  the  mechan¬ 
ism,  particularly  in  the  combustion  environment,  made  estimates  highly 
uncertain. 

Particulate  damping  was  found  to  be  a  strong  dissipative  factor,  particu¬ 
larly  at  high  frequency.  Droplet  deformation  was  completely  insignificant. 
The  contri nut  ion  due  to  surface  modes  of  vibration  was  larger  than  that 
for  droplet  deformation,  but  in  absolute  magnitude  was  also  insignificant. 

In  the  absence  of  experimental  data  concerning  the  shattering  of  drops 
near  the  burning  surface,  ns  well  as  a  model  for  their  subsequent  behavior, 
droplet  shattering  as  a  dissipative  mechanism  could  only  be  bracketed 
between  wide  limits,  with  the  upper  limit  of  significant  relative  magnitude. 
The  effects  of  aluminum  on  combustion  interaction  of  the  burning  surface 
with  the  acoustic  field  was  found  to  he  a  completely  intractable  problem 
within  the  scope  of  the  present  program;  in  consequence ,  no  quantitative 
estimate  of  iis  associated  energy  dissipation  rate  could  be  made. 

In  addition,  tin  effects  of  aluminum  on  nonacoustic  instability  were 
considered.  An  investigation  of  nonacoustic  instability  produced  several 
results  of  considerable  importance,  although  not  necessarily  v'ith  respect 
to  the  role  of  aluminum.  Analysis  showed  that  the  nonacoustic  problem 
could  he  developed  using  conventional  acoustic  methods  with  the  nonacoustic 
mode  being  a  zeroth-longitudi nal  acoustic  mode  for  the  cavity.  With  this 
basis,  an  acoustic  instability  combustion  model  was  applied  to  the  descrip¬ 
tion  of  the  nonacoustic  instability  data  reported  by  several  investigators. 
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Highly  .satisfactory  agreement  was  obtained.  Finally,  the  information 
developed  was  applied  to  the  effects  of  aluminum  in  the  nonacoustic 
frequency  range  to  propose  a  source  of  the  periodicity  in  aluminum 
shedding  observed  in  experimental  investigations. 

The  results  of  the  analytical  task  made  the  need  for  experimental  measure¬ 
ments  evident.  A  method  was  needed  to  measure  the  combustion  response 
(response  function)  of  the  burning  surface  for  a  highly  aluminized  propell¬ 
ant.  Secondly,  any  changes  in  the  response  function  needed  to  be  compared 
to  the  particulate  damping  for  that  propellant,  to  assess  the  importance  of 
these  changes  in  the  response  function.  Therefore,  measurements  of  the 
particulate  damping  were  also  needed.  These  measurements  were  made 
employing  a  T-hurner  in  which  a  pressure  pulse  could  be  generated. 

Although  this  method  was  not  fully  satisfactory  with  respect  to  the 
accuracy  obtainable,  there  is  at  present  no  other  method  capable  of 
producing  the  necessary  information.  Using  this  technique,  measurements 
were  made  on  propellants  containing  up  to  15-percent  aluminum  at  mean- 
pressure  levels  to  600  psi  and  over  a  frequency  range  of  200  to  2000  ops. 


Satisfactory  results  were  obtained.  1+  was  found  that  particulate  damping 
was  the  dominant  stabilizing  mechanism.  Measurable  particulate  damping 
was  obtained  even  at  the  lowest  frequencies.  The  response  function  did 
not  appear  to  be  greatly  altered  by  the  presence  of  aluminum;  however, 
the  data  consistently  indicated  that  the  aluminized  burning  surface  is  more 
unstable  than  the  corresponding  unaluminized  surface.  In  spite  of  limita¬ 
tions  of  the  experimental  method,  the  consistency  of  the  data  suggests 
that  the  conditions  exist.  These  results  are  of  considerable  importance 
because  they  represent  tne  only  data  available  for  highly  aluminized 
propellants . 
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ANALYSIS 

i' r  of  ;  hei  »"-r»  nn  Day  bo  cited  as  possible  mechanisms  by  which  inclusion 
n»i  .•;«?,  powder  in  a  solid  propellant  formulation  introduces  energy 
s  or  sinks  during  unstable  combustion.  For  some  of  these  (e.g. , 
part  i  eutntc  damping;  indications  existed  before  initiation  of  the  reported 
•te’y  that  their  role  vas  important.  In  most  cases,  however,  no  evidence 

Iti  x!n>  first  phase  of  the  investigation,  the  following  possible  mechanisms 

of  !:*.r;p :  nr  were  eon  si  tiered  analytically: 

1 ,  Pu r t  i  c  u  1  n  t  e  d  amp  ing 

2.  Deformation  of  droplets  at  the  burning  surface 

Droplet  surface  modes 

i,  Droplet  shattering 
>.  Combust  i. on  interaction 

For  M'h  of  those  it  was  attempted  to  develop  a  quantitative  expression 
for  the  change  of  energy-vi th-time  associated  with  the  mechanism.  The 
occurrence  of  some  of  these  phenomena  (c.g. ,  droplet  shattering)  in  the 
vi c : n i ty  of  the  burning  surface  considerably  complicates  the  analysis. 

If  th<>  latter  is  to  hr  kept  within  the  bounds  of  traetnbility ,  simplifica- 
ti  ■:!;?  arc  r r’cessary  which  sometimes  lead  to  results  which  are  valid  only 
■s  eider  ("■  f  rr  i  Piido  calculations. 

fyh  .  Muer !y ,  v:t  y-kt .in,-;  factors  were  calculated  for  each  mechanism  so  that 
<•  i  v-  :  rrr;!,)';-.: ‘y-,0  6  h*>  made.  This  required  specific  assumptions  to  be 
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made  concerning ,  in  particular,  the  ratio  of  the  number  of  aluminum  droplets 
on  the  burning  propellant  surface  to  those  in  the  produc t-gas  volume.  Tt 
was  then  possible  to  directly  estimate  the  relative  influence  of  the 
dissipation  factors  considered.  In  spite  of  the  approximations  made,  in 
most  cases  contributions  differ  by  a  great  enough  amount  to  permit  clear 
resolution  of  this  relative  influence.  Where  this  resolution  has  not  been 
possible  by  the  analytical  techniques  described  in  this  section,  further 
evidence  on  d.imping  was  obtained  from  the  experimental  phase  of  the  program. 


PARTICEUTE  DAMPING 

Particulate  clamping  is  indicated  as  a  major  effect  in  the  role  played  by 
powdered  aluminum  in  suppressing  acoustic  instability.  Consequently,  it 
is  convenient  to  discuss  this  process  first  and  use  it  as  a  basis  for 
comparison  with  the  contributions  of  other  damping  processes.  For  these 
reasons,  as  well  as  its  greater  tractability  than  many  of  the  mechanisms 
treated  subsequently ,  particulate  damping  has  been  analyzed  in  considerable 
detai 1 . 

The  theory  of  particulate  damping  has  been  analyzed  by  Sewell  (Ref . 9  ), 

Lamb  (itef.  10  ) ,  Epstein  and  Carhart  (Ref.  1 1) ,  Blair  (Ref.  12),  Chow  (Ref.  13 
Soo  (Ref.  14  ),  and  most  recently  by  Temkin  (Ref.  13  ).  The  results  of 
these  investigations  appear  to  be  essentially  in  agreement  with  the 
exception  of  Soo  (Ref.  14  ).  Soo  finds  a  small  difference  (about  a  factor 
of  two)  between  his  results  and  those  of  Epstein-Carhart  (Ref.  ll)  for  a 
particular  numerical  example.  However,  Soo's  model  is  equivalent  to  that 
of  Temkin  (Ref.  13)  over  the  range  of  interest  but,  contrary  to  Soo, 

Temkin  finds  excellent  numerical  agreement  between  his  results  and  the 
Epstein-Carhart  results.  In  the  current  investigation,  the  Temkin  analysis 

( 
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was  verified  and  excellent  numerical  agreement  was  also  found  with  the 
Epstein-Carhart  results.  Therefore,  the  disparity  was  discounted  as 
unimportant.  In  the  current  program  the  Epstein-Carhart  equations  were 
used  for  most  calculations.  The  Temkin  equation  was  used  for  comparative 
purposes  and  for  some  calculations. 
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The  theory  has  had  substantial  experimental  verification,  although  the 
data  are  limited.  The  data  of  Zink  and  Delsasso  (Ref.  32  )  agree  very 
well  with  Epstein-Carhart  theory,  over  the  range  in  which  data  were  ob¬ 
tained.  The  recent  data  of  Temkin  and  Dobbins  (Ref.  33)  agree  very  well 
with  the  Temkin  treatment  (and,  presumably,  Epstein-Carhart  theory  as 
well)  over  a  broad  range  of  conditions.  Earlier  data  by  Dobbins  and 
Temkin  (Ref.  34  )  also  agree  with  Epstein-Carhart  theory.  Therefore,  it 
appears  that  particulate  damping  is  rather  well  understood. 

The  most  well-known  analysis  of  particulate  damping  is  the  theory  of 
Epstein  and  Carhart  (Ref.  11  ).  They  have  solved  the  wave-scattering 
problem  for  a  plane  wave  impinging  upon  a  single  droplet,  including 
viscous  and  thermal  effects.  The  time-average  rate  of  energy  dissipation 
was  then  calculated  from  the  rate  of  entropy  generation  in  the  entire 
volume, both  within  and  without  the  sphere.  The  assumption  was  made  that 
the  rate  of  dissipation  for  many  droplets  may  be  found  by  adding  the 
dissipation  rate  for  each  droplet.  Thus  they  calculate  an  attenuation 
coefficient  for  an  aerosol.  The  reduced  attenuation  coefficient  is 
given  by  (Ref.  34  ): 


a  a 
o 


taS  =  6  =  6  {  4  (1  +  y)  F  (y)  +  ^  (1  +  Prl/2  y)  G  (y)}  (l) 


m  * 
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where 

F  (y) 
G  (y) 


(2) 

(3) 


The  nomenclature  is  shown  in  Appendix  A. 

The  approach  taken  by  Temkin  (Ref.  15  )  is  rather  different  from  that 
used  by  Epstein  and  Carhart.  He  writes  the  conservation  equations  for 
a  particle-laden  pas.  The  assumptions  were  made  that  Stoke's  law  may  be 
used  to  describe  the  force  on  a  particle  and,  secondly,  that  the  heat 
transport  to  a  particle  may  be  described  by  the  expression  for  heat  con¬ 
duction  to  a  sphere  in  a  stagnant  fluid.  The  linearized  conservation 
equations  were  solved  for  one-dimensional  behavior.  For  very  small  mass 
concentrations  of  particles,  he  finds: 


a  = 


CO  T 


d 


1  +  60  r 

a 


(4) 


If  the  thermal  dissipation  is  negligible  relative  to  the  viscous  dissipa¬ 
tion,  then  for  finite  mass  concentrations 


Cm  l 

2(\  +  o?  rd2)J 

Detailed  calculations  were  made  with  both  the  Epstein-Carhart  equation 
and  the  Temkin  equation  as  noted  above.  A  numerical  example  for  low 
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concentrations  ae  calculated  by  each  equation  is  shown  in  Fig.  1  .  The 

calculations  were  made  for  Al^O^  droplets  suspended  in  combustion  gases. 

The  numerical  agreement  is  very  good  except  at  large  values  of  y,  i.e., 

,IC  r2  1/2 

(7 7”p~)  *  Perhaps  the  disparity  in  this  region  can  be  attributed  to 

deviations  from  the  simple  Stokes-lav  drag  relationship  expected  in  the 
region,  fortunately ,  the  Temkin  equation  is  adequate  for  the  combustion- 
instability  case  since  the  particle  sizes  encountered  are  small. 

In  Fig.  2  numerical  comparison  is  shown  between  the  viscous  contribution 

and  tho  total  reduced  damping  coefficient  as  calculated  by  the  Temkin 

equation.  It  is  apparent  that  the  thermal  contribution  is  Bmall  for  this 

case.  Experimentally,  in  combustion  instability  studies,  the  temporal 

attenuation  coefficient  is  the  quantity  usually  measured.  The  Epstein- 

Carhart  curve  shown  in  Fig.  1  has  been  converted  to  such  a  coefficient 

and  the  results  are  plotted  in  Fig.  3  for  several  particle  sizes  and 

a  mass  ratio,  C  ,  of  one-tenth.  In  these  terms  the  effectiveness  of 
m 

particulate  damping  is  very  apparent. 

It  is  useful,  in  the  discussion  of  experimental  results,  to  consider  the 
mean  particle  size  to  be  used  in  these  equations  if  a  distribution  of 
particle  sizes  is  present.  The  mass  ratio,  C^,  may  be  estimated  from 
the  propellant  composition  and,  therefore,  is  preoumed  to  be  known. 

The  thermal  contribution  is  assumed  to  be  negligible.  It  is  simplest 
to  deal  with  the  Temkin  equation;  the  Epstein-Carkart  equation  is  dealt 
with  in  an  identical  manner.  Integrating  over  the  size  distribution, 
it  is  found  that 
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Figure  1.  Comparison  of  Fpstein-Carhart  and  Tpmkin  Results 
for  Particulate  Damping 
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Figure  3 


Particulate  Damping  in  Terms  of  the  Temporal 
Damping  Constant  According  to  Epstein- 
Carhart  Theory 
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Two  simple  cases  occur.  If  the  distribution  is  such  that  the  integrand 
is  negligible  foi 
k  r  <<c  1 ,  then 


is  negligible  for  large  sizes,  i.c.,  the  integrand  is  non-negligible  for 
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i,  j  =  0,  1,  2,  3, 


(10) 


it  is  apparent  that  the  D  mean  diameter  is  appropriate  for  small 

53 

particles. 

2  4^ 

If  the  integrand  is  non-negligible  for  r  >:>  1 ,  then  the  D.^  mean 
diameter  is  found  for  large  particles.  From  the  Epstein-Carhart  equation 
it  is  also  found  that  should  be  used  for  small  particles,  D3i  for 
large  particles,  and  D,„  for  very  large  particles.  These  considerations 

J*- 

were  used  in  examining  the  experimental  results  discussed  below. 
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The  rate  of  energy  dissipation  due  to  particulate  damping  may  be  obtained 
rather  directly.  Assuming  the  thermal  contribution  to  be  negligible,  it 
was  found  from  the  Epstein-Carhart  theory  that,  for  each  particle, 


Epart  *  3  n  T  v  (i  ♦  y)  y F  (y) 


for  small  particles,  this  reduces  to 


(11) 


4  IT 


r5  t? 


part 


27  *  62  P  aQ2 


a2 

P 


(12) 


where  $  is  the  amplitude  of  the  incident  pressure  wave.  Equation  12 
will  be  used  to  make  comparisons  for  specific  cases  after  expressions 
have  been  described  for  other  dissipative  processes. 

ENERGY  DISSIPATION  BY  DEFORMATION  OF 
DROPLETS  ON  THE  BURNING  SURFACE 

Perhaps  the  simplest  technique  available  fcr  calculation  of  the  rate  of 
degradation  of  acoustic  energy  is  described  by  Landau  and  Lifshitz 
(Ref.  17).  From  consideration  of  the  rate  of  entropy  generation,  they 
obtain  the  following  expression; 


e =  -  <  1  r  J  (7T)2  dv  +  **  J  % dv} 


(13) 


0  V 


where  the  integration  is  carried  out  over  the  entire  volume  and  <  >  in¬ 
dicates  a  time  average.  The  convenience  of  Eq.  13  is  that  E  to  a  first 

f 

approximation,  can  be  calculated  from  the  solutions  to  the  acoustic  problem 
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with  losses  neglected,  since  the  losses  contribute  only  second-ordpr 
effects  in  general.  This  approximation  was  used  to  calculate  the  dis¬ 
sipation  in  the  various  oscillatory  modes  of  a  droplet.  A  similar  approx¬ 
imation  was  used  by  Cantrell  (Ref.  18  ). 

Two  cases  were  considered:  (l)  a  drop  located  near  a  velocity  node  and 
(2)  a  drop  exposed  to  a  transverse  acoustic  velocity.  In  each  case  the 
drop  was  approximated  as  spherical.  In  the  first  case,  only  radial  be¬ 
havior  need  be  considered.  No  problem  is  encountered  in  attaining  a 
solution  to  the  wave  equation  for  this  case.  The  viscous  contribution 
to  the  dissipation  was  found  to  be 

E  .  =  16  IT  JU.D 

vis  l 

v  ,  where  terms  smaller  than 

tribution  was  found  to  be 

\  A2  .  _4 

Etherm.  =  4  T* J  *Vis  "  10  *Vis 

'  * 

which  is  negligible.  Eq.  14  will  be  used  later  for  specific  caces. 

For  the  droplet  exposed  to  a  transverse  acoustic  velocity,  some  difficulty 
was  found  in  handling  the  velocity  distribution  in  the  boundary  layer. 

This  distribution  was  finally  handled  in  an  approximation.  A  solution 
was  obtained  in  a  manner  similar  to  that  used  by  Rayleigh  (Ref.  19  ). 

The  solution  was  found  to  exhibit  predominantly  radial  motion,  the 
nonradial  portion  of  the  solution  being  on  the  order  of  10~^  times  the 
radial  portion.  This  result  may  be  due  to  the  imposed  approximations; 
however,  on  the  basis  of  the  approximate  computation,  it  was  concluded 
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that  the  contribution  to  the  dissipation  was  of  the  same  order  as  that 
found  for  a  drop  located  near  a  velocity  node. 

ENERGY  DISSIPATION  IN  THE  SURFACE 
MODES  OF  A  DROPLET 

The  resonant  frequencies  for  aluminum  droplets  oscillating  in  their 
lowest  surface  modes  are  in  the  range  of  instability  frequencies  of  in¬ 
terest,  for  not  unreasonably  large  droplets.  This  suggests  that  a  high 
dissipation  may  be  obtained  from  those  droplets.  The  dissipation  rate 
was  again  calculated  employing  Eq.  13  . 

The  solution  to  the  surface-mode  problem  is  given  by  Rayleigh  (Ref.  19), 
and  Landau  and  Lifshitz  (Ref.  17  ),  for  example.  By  integration  of 
Eq.  13  with  this  solution  it  was  found  that,  for  the  lowest  mode: 

1  2 
31  IT  D  J  A 

*  -  — sr —  he) 

The  lowest  resonant  frequency  is  given  by 
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DISSIPATION  BY  PERIODIC  DROPLET  SHATTERING 


If  droplets  are  shattered  in  a  periodic  manner,  acoustic  energy  can  be 
converted  to  surface  energy  of  the  drops.  The  increase  in  energy  is 
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proportional  to  the  increase  in  surface  area.  Assuming  a  droplet  to  be 
broken  into  several  equally  sized  droplets,  the  energy  increase  is 

A  E  =  CT  77  Dg2  (n2/3  -  l)  (18) 

where  n  is  the  number  of  droplets  formed  from  the  initial  droplet.  The 
difficulty  that  arises  in  this  case  is  in  estimating  the  frequency  of 
shattering. 

Various  studies  of  droplet  shattering  have  been  made  (Ref. 35  through 37)« 

In  these  studies  the  interest  was  the  behavior  of  drops  exposed  to  open 
velocity  fields,  e.g.,  the  flow  field  following  a  shock  wave.  Unfortunately, 
this  is  a  different  situation  from  that  encountered  by  a  droplc:  located 
on  a  surface.  These  studies  have  shown,  however,  that  the  principal  source 
of  breakup  is  associated  with  the  relative  velocity  between  the  gas  and 
droplet,  rather  than  impulsive-pressure  forces.  This  indicates  that 
shattering  of  a  drop  on  a  surface  will  not  be  important  unless  it  is 
exposed  to  a  tangential  velocity.  It  is  likely  that  the  drops  will  be 
swept  from  the  surface  before  shattering  occurs  under  these  circumstances. 

The  drop  removal  will  involve  an  energy  exchange,  in  the  form  of  surface 
energy  and  kinetic  energy,  but  this  is  a  considerably  less  dissipative 
effect  than  shattering  a  droplet  into  several  "sub-droplets."  Because 
shattering  studies  indicate  no  means  of  estimating  the  number  of  drops 
formed  from  a  single  drop  or  the  frequency  of  drop  shattering,  the  com¬ 
parison  in  a  subsequent  section  is  made  with  assumed  values  for  these 
parameters. 
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EFFECTS  OF  ALUMINUM  ON  COMBUSTION  INTERACTION 
OF  THE  BURNING  SURFACE 

It  seems  possible  that  agglomerates  of  aluminum  on  the  burning  surface 
could  physically  block  interaction  with  the  acoustic  field.  However, 
the  time-average  buring  rate  of  solid  propellants  is  not  greatly  affected 
by  the  addition  of  aluminum.  This  suggests  that,  on  the  average,  the 
non-aluminum  portion  of  the  propellant  burns  about  the  same  as  a  non- 
aluminized  propellant.  From  this  it  appears  unlikely  that  one  can  even 
crudely  estimate  the  effect  of  aluminum  in  this  case.  Such  estimates 
must  await  analysis  of  the  heterogenous  combustion  problem  involving  the 
aluminum.  This  is  an  exceedingly  complex  problem,  and  appears  totally 
intractable. 

In  the  same  category  fall  those  effects  caused  by  the  aluminum  which  in¬ 
volve  disturbance  of  the  energy  feedback  to  the  surface.  If  the  role 
played  by  these  effects  in  combustion  instability  is  to  be  determined  it 
must  be  done  through  careful  experimentation,  allowing  delineation  of 
their  influence.  The  measurement  of  surface  and  chamber-volume  effects 
which  constituted  the  experimental  part  of  this  program,  allows  an  esti¬ 
mate  of  these  influences  to  be  made. 


Comparison  of  Damping  Mechanisms 

To  compare  particulate  damping,  the  principal  volumetric  effect  associated 
with  the  presence  of  the  aluminum,  with  aluminum  effects  on  or  near  the 
burning  surface,  it  is  necessary  to  weight  each  contribution  by  the 
number  of  droplets  of  each  type.  Therefore,  the  number  of  droplets  in 
the  main  volume  relative  to  the  number  on  the  surface  must  be  estimated. 
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Photographs  of  the  burning  surface  of  heavily  aluminized  propellants 
have  shown  the  surface  to  be  nearly  covered  with  molten  aluminum,  at 
least  at  times.  To  approximate  this  situation,  the  surface  was  assumed 
to  be  covered  by  spheres  which  touched  each  other  and  whose  line  of 
centers  formed  a  square  pattern.  The  number  of  drops  on  such  a  surface 
is  given  by 


N 

s 


(19) 


The  number  of  droplets  suspended  in  the  gas  phase  can  be  obtained  from 
the  mass  ratio,  C^,  assuming  its  value  to  be  known.  From  the  definition 
of  C^,  the  number  of  droplets  in  found  to  be 

6  C  6  V 

N  =  - 5-=-  (20) 

V  TTD  3 
v 


The  value  of  C  was  calculated  by  assuming  all  of  the  aluminum  to  be 
oxidizer  to  Al^O^. 


For  a  motor  with  a  cylindrical  perforation  and  burning  radially,  the 
ratio  of  numbers  of  drops  is  found  to  be 


N  ,  D  2  D  C 
v  3  _s _ 
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6 

m 


(21) 


For  a  port  diameter  of  six  inches,  an  aluminum  loading  of  15  weight 
percent  (or  C  =  0.396),  and  6  =  5-3  (l0  ),  the  values  shown  in  Table  1 

were  calculated. 
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TABLE  1 


NUMBER  OF  DROPS  IN  GAS  VOLUME  RELATIVE 
TO  NUMBER  OF  DROPS  ON  BURNING  SURFACE 


D  ,  microns 

V 

D  ,  microns 
s 

NA 

0.1 

100 

1.2  (108) 

1.0 

100 

1.2  (105) 

10.0 

100 

1.2  (102) 

0.1 

1000 

O 

H 

O 

r— 1 

CM 

1.0 

1000 

1.2  (107)  - 

10.0 

1000 

1.2  (104) 

With  the  values  in  Table  1  ,  the  appropriately  weighted  contribution 
from  particulate  damping  was  computed.  Choosing  a  pressure  amplitude 
of  10  psi  and  a  frequency  of  1000  cps  as  representative,  the  values 
shown  ir\  Tuble  2  were  calculated  using  Eq.  14  . 


TABLE  2 

WEIGHTED  RATE  OF  ENERGY  DISSIPATION 
DUE  TO  PARTICULATE  DAMPING 


D  ,  micron 

V 

D  ,  micron 

s 

<W  Eparf 
(ergs/drop-sec ) 

0.1 

100 

1.4  (101) 

0.1 

1000 

1.4  (103) 

1.0 

100 

1.4  (103) 

1.0 

1000 

1.4  (105) 

10.0 

100 

1.4  (i03) 

10.0 

10U0 

1.4  (107) 

on 
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This  table  provides  the  reference  values  for  other  estimates.  No  attempt 
has  been  made  to  cover  a  complete  range  of  conditions. 

The  dissipation  rate  associated  with  droplet  deformation  was  estimated 
by  use  of  Eq.  14  .  With  the  following  properties: 

nl  cl  -  ££  =  7,5  (1q11)  dyne/cm2 
c^  ~  5*5  (106)  cin/sec 
JU(,  ~  2  (lC  2)  poise 
f  =  lUOO  cps 
"p"  =  10  psi 

the  results  shown  in  Table  3  were  obtained. 

TABLE  3 

RATE  OF  ENERGY  DISSIPATION  DUE  TO 
DROPLET  DEFORMATION 

D^ ,  microns  E  (erg/drop-sec) 

100  1.2  (10~5) 

1000  1.2  (10~4) 

Therefore,  droplet  deformation  appears  to  be  insignificant  relative  to 
particulate  damping,  i.e.,  Table  2  . 
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The  contribution  due  to  surface  modes  of  oscillation  was  calculated  from 
Eq.  lo  .  The  assumptions  made  in  calculating  ft  for  this  case  were:  (l) 
that  resonance  is  attained  by  the  droplet,  (2)  that  the  maximum  disturbance 
of  the  radius  of  the  drop  is  10  percent  of  its  mean  value,  and  (3)  the 
surface  tension  is  840  dyne/cm  (appropriate  for  aluminum  in  contact  with 
air  at  700  C).  At  1000  cpe,  it  was  found  that  ft  =  7*3  erg/drop-sec.  The 
corresponding  drop  diameter  was  830  microns. 

Thus,  at  the  chosen  conditions,  the  surface  modes  are  not  comparable  to 
particulate  damping.  They  appear  to  be  insifnificant ,  in  spite  of  their 
relatively  greater  magnitude  with  respect  to  droplet  deformation. 

It  is  more  difficult  to  make  meaningful  estimates  in  the  case  of  droplet 
shattering.  One  can  only  assume  values  for  the  frequency  of  shattering 
and  the  number  of  new  drops  formed  from  an  initial  drop.  If  a  single, 
1000-micron  drop  is  shattered  into  twenty  "sub-droplets,"  Eq.  18  predicts 
175  ergs  of  energy  will  be  absorbed.  The  importance  of  shattering  as  a 
dissipation  phenomenon  then  becomes  dependent  upon  the  shattering  frequency, 
since  the  energy  associated  with  the  individual  shattering  event  is  sub¬ 
stantial.  As  already  noted,  there  is  no  exnerimental  evidence  to  draw 
upon.  Thus,  shattering  r '  the  surface  may  not  occur  at  all.  It  is  likely 
that,  if  it  does  occur,  oiuy  a  small  number  of  drops,  on  the  average,  would 
shatter  during  each  acoustic  oscillation  (this  number  may  be  fractional). 

The  dissipation  rate  for  the  shattering  process  is  the  product  of  the 
energy  loss  per  shattering  event,  the  mean  number  of  events  per  acoustic 
oscillation,  and  the  acoustic  frequency.  It  is  clear  from  the  discussion 
that  the  lower  limit  is  negligibly  small.  Under  optimum  conditions  for 
this  process,  however,  it  may  compare  favorably  with  particulate  damping 


no 
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as  a  dissipation  factor,  (it  must  be  borne  in  mind  that  other  factors, 
not  at  all  understood  at  present,  also  enter  into  overall  consideration 
of  the  droplet  shattering  mechanism.  For  example,  what  further  role,  if 
any,  is  played  by  droplets  resulting  from  the  shattering  of  a  surface 
drop.  Also,  after  shattering  occurs  does  reformation  through  coalescence 
occur  to  a  significant  extent). 

The  analysis  presented  indicates  that  particulate  damping  is,  in  general, 
the  dominant  dissipation  mechanism  in  the  high  frequency  range,  and  for 
the  processes  considered.  It  is  possible,  however,  that  the  effects  of 
the  aluminum  on  combustion  interaction  with  the  acoustic  field  are  im¬ 
portant  to  stability.  A  similar  statement  may  be  made  concerning  droplet 
shattering.  These  dissipation  mechanisms  must  be  accepted  as  insufficiently 
accessible  to  analysis  of  the  desired  precision,  so  that  further  success 
with  this  problem  is  possible  only  through  experiment.  The  results  of  the 
experimental  portion  of  this  program  suggest  that  these  surface  effects 
may  indeed  play  a  role  in  characterizing  the  "stability"  of  a  propellant, 
although  a  role  subordinate  to  that  of  particulate  damping. 

Before  considering  the  experimental  phase  of  the  program,  a  brief  digression 
will  be  made  to  consider  a  portion  of  the  analysis  which  was  devoted  to 
nonacoustic  instability. 


THE  EFFECTS  OF  ALUMINUM  ON  NONACOUSTIC  INSTABILITY 

Nonacoustic  instability  has  been  studied  rather  extensively  (Ref.  1, 

20,  22,  and  28).  These  studies  suggest  that  the  inclusion  of  powdered 
aluminum  may  be  a  destabilizing  influence,  i.e.,  the  aluminum  has  a  det¬ 
rimental  effect  on  combustion  instability.  In  support  of  this  suggestion 
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is  the  observation  that  aluminum  is  shed  quasi-periodically  from  the 
surface.,  with  a  frequency  in  the  range  of  no.oacouatic  instability,  when 
burned  in  the  open  or  in  a  closed  bomb.  Therefore,  it  is  of  ir+erest 
to  study  the  effects  of  aluminum  on  nonacouBtic  as  well  as  acoustic 
instability. 

To  learn  something  about  these  effects,  an  analytical  investigation  of 
nonacoustic  instability  was  carried  out.  The  analysis  is  described  in 
Appendix  B.  Several  important  results  were  obtained,  although  most  of 
these  are  not  directly  applicable  to  the  effects  of  aluminum.  It  was 
found  that  the  nonacoustic  problem  could  be  developed  using  conventional 
acoustic  methods,  the  nonacoustic  mode  being  a  zeroth-longitudinal  acoustic 
mode  for  the  cavity.  This  suggested  the  application  of  an  acoustic- 
instability  combustion  model  to  the  nonacoustic  case.  The  Denison-Baum 
(Ref.  25)  model  was  found  to  describe  the  nonacoustic  data  of  Beckstead 
(Ref.  22  )  very  well.  In  addition,  a  detailed  comparison  between  the 
Denison-Baum  model  and  the  acoustic  (T-burner)  instability  data  of  Horton 
(Ref. 27  )  was  “ade.  Again,  rather  good  agreement  was  attained.  No  pre¬ 
vious  comparisons  of  this  type  have  been  reported.  The  combined  model 
was  found  to  explain  the  "preferred  frequency"  reported  by  NOTS  (Ref.  28  and 
29)  and  phase-angle  data  between  the  oscillatory  pressure  and  light 
emission,  also  reported  by  NOTS  (Ref.  29). 

With  reference  to  the  effects  of  aluminum,  these  results  suggest  two 
things.  First,  since  the  response  function  (Appendix  B)  exhibits  a 
maximum  in  the  nonacoustic  frequency  range,  it  appears  at  least  possible 
that  the  quasi-periodic  shedding  of  aluminum  coincides  with  this  maximum. 
This  possibility  is  arrived  at  by  the  following  argument:  Pressure  dis¬ 
turbances  generated  by  the  combustion  will  be  reflected  back  toward  the 
burning  surface  by  the  dense  cloud  of  aluminum  leaving  the  surface. 
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*'  tre,,,“Cy  '""'•“'"‘J  “f  «■«  pressure  disturbance  which 

arc  near  the  »s,iM  in  f  win  be  amplified  to  .  greater  extent  than 

o  hers.  Thus,  lacking  an  externally  imposed  frequency,  the  combustion 

8  °  ten<l  t<’Varl1  ltB  resP°nse,  and,  presumably,  so  will  the 

shedding. 

The  second  conclusion  of  this  vorh  is  that,  if  the  addition  of  aluminum 

has  .  dxstaMlizing  in  W-,  then  the  influence  he  reflected  in 

e  e  (€)-  That  ia,  Re  (  )  for  an  aluminized  propellant  must  be  greater 
m  the  low  frequency  Tang,  than  that  for  greater, 

g  ’  that  for  th<-  corresponding  nonaluminized 

propellant. 


Alth  h  en'3  °f  thl'  Pr°g‘“’  are  °f  i”  ‘his  sense. 

houEh  he  experimental  data  are  not  sufficiently  precise  to  warrant 

irm  conclusions,  the  data  consistently  indicate  that  the  Re  («)  i,  in 

creased  by  the  addition  of  aluminum.  In  the  high-freqnency  rjige,  the 

increase  in  volumetric  damping  overwhelm,  the  increase  in  surface  amplifi- 

catxon.  In  the  nonacoustic  fremipnrv  . 

tic  irequency  range,  however,  particulate  damp¬ 
ing  is  ineffectual  and  the  am*.--— *■  „  .  _  ,U\ 

-P  v  iucrease  in  Re  (f)  is  very  important. 
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EXPEEUMENTAL  STUDY 


it  i  apparent  tha  t  an  experimental  method  was  needed  to  measure  the  ef¬ 
fect:  of  aluminum  on  combustion  instability.  Accurate  data  are  needed 
on  the  volumetric  ivgas  phase)  damping  and  burning-surface  response  function 
for  highly  a  its:  faired  propellants  and  the  corresponding  non— aluminized 
propellants.  Unfortunately,  such  data  has  in  the  past  been  inaccessable 
because  of  the  difficulty  in  making  the  appropriate  measurements  and  de- 
v 5 sing  the  necessary  experiments.  The  usual  T-burner  techniques,  which 
ore  fairly  v ideiy  t**ed  to  study  acoustic  instability  of  nonmetalized 
propellants,  are  not  applicable  because  of  the  high  stability  of  the 
aestniized  system.  Horton  (Bef.  6)_  used  a  conventional  T-burner  to  study 
the  effects  of  aluminum  but  he  was  unable  to  employ  aluminum  loadings 
greeter  than  1.5  freight  percent.  Although  his  study  was  productive ,  the  j 

airrinuS;  loadings  he  employed  were  an  order  of  magnitude  too  low  to  be  ; 

applicable  to  conventional  propellants.  A  second  technique  is  the  shock- 
tihi  metnod  (ft-f.  JO)  which  has  been  utilized  at  fiocketdyne.  However, 
lit i a  technique  provides  no  information  on  particulate  damping,  apparently 
the  most  important  damping  process.  .Also,  experimental  and  computational  ‘ 

difficulties,  prevent  easy  application  of  the  method.  The  technique  which  ; 

.tt'f'Tnod  applicable  to  obtaining  the  requisite  information  was  a  pulsed 

i 

.has  shown  (Kef.  JT)  that,  if  a  pressure  pulse  is  properly  intro-  I 

dccvl.  a  T-bumer  ecu  be  used  to  study  propellants  too  stable  to  drive  5 

oscillations.  The  pressure  pulse  must  be  such  that  transient  pressure  | 

i  llatioas  are  generated  in  the  T-burner.  ?■ 


I 


preceding  page  blank 
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Further,  measurements  of  the  rate  of  decay  of  these  oscillations  provide 
sufficient  information  to  calculate  the  response  function  of  the  burning 
surface  and  the  volumetric  damping.  Even  though  this  method  suffers  from 
certain  limitations  it  was  applied  as  the  best  method  currently  available. 

The  usual  T-buraer  technique  is  to  measure  the  growth  constant  from  the 
rate  of  oscillatory  development  at  the  beginning  of  a  test  and  the  decay 
constant  from  the  oscillatory  decay  following  burnout  (Eq.  Ref.  38,  39 j 
and  40).  These  two  kinds  of  derived  data  may  be  interpreted  to  calculate 
the  response  function  for  the  burning  surface.  The  decay  constant  is 
taken  to  represent  the  overall  damping  in  the  system.  In  a  pulsed  T- 
burner,  the  oscillations  do  not  develop  because  the  damping  exceeds  the 
gains.  Therefore,  after  the  burning  is  well,  established,  a  pulse  is 
introduced  which  excites  oscillations  in  the  cavity.  These  oscillations 
decay  because  of  the  excess  damping.  However,  this  rate  of  decay  may  be 
interpreted  in  a  manner  identical  with  the  usual  growth  constant.  A 
second  pulse  is  introduced  immediately  after  burnout  to  measure  the  over¬ 
all  damping.  This  damping  will  include  whatever  damping  is  introduced 
by  the  pulsing  device.  Thus,  the  interpretation  of  pulsed  T-burner  data 
is  carried  out  in  a  manner  essentially  identical  with  conventional  T- 
burner  data. 


EXPERIMENTAL  SYSTEM 

The  T-burner  used  in  this  study  was  similar  to  that  used  by  other  investiga¬ 
tors.  A  sketch  of  the  motor  is  shown  in  Fig.  4  ,  The  motor  was  made  of 
steel  tubing  with  dimensions  of  1.5-inch  ID  by  0.5-inch  wall  thickness 
(nominally).  Interchangeable  sections  were  made  si  that  motor  lengths 
ranging  from  41  inches  to  9  inches  (without  propellant)  were  available. 

These  lengths  provided  instability  frequencies  of  200  cps  to  somewhat  more 
than  2000  cps. 
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Squibs  were  used  to  generate  a  pressure  pulse  when  necessary.  Two  squibs 
were  mounted  in  the  curved  wall  of  the  motor.  For  proper  pulsing  it  was 
necessary  to  have  a  positive  indication  of  propellant  burnout  in  the  motor. 
A  thermocouple  was  mounted  under  one  propellant  sample  to  serve  this 
function. 

A  Kistler  pressure  transducer  (model  603A)  was  used  to  measure  the  oscil¬ 
latory  pressure.  Except  for  the  shortest  motor,  the  transducer  was  mounted 
in  the  curved  wall  of  the  motor,  about  2-3/4  inches  from  one  end.  In  the 
shortest  motor,  the  transducer  was  mounted  under  the  propellant  sample 
in  one  end  closure. 

Propellant  samples  were  bonded  into  each  end  of  the  motor  with  an  epoxy 
adhesive.  The  propellant  was  cast  into  phenolic  sleeves  (l/l6-inch  wall 
thickness)  and  then  cut  on  a  lathe  into  l/2-inch  lengths.  The  phenolic 
sleeves  facilitated  bonding  and  prevented  burning  on  the  sides  of  the 
sample,  which  can  significantly  alter  results.  Lathe  cutting  provided 
a  uniform  and  reproducible  surface,  as  well  as  an  accurate  length. 

Ignition  was  obtained  by  burning  a  small  amount  of  propellant  shavings 
near  the  propellant  surface.  The  shavings  were  ignited  by  burning  a 
1-1/2-inch  (nominal)  length  of  exploding  bridge  wire  (Pyrofuze).  The 
bridge  wire  was  soldered  to  the  electrical  leads  of  an  expended  squib 
(initiator)  which  eliminated  sealing  problems  and  facilitated  replacement. 

The  motor  was  close-coupled  to  a  50-gallon  surge  tank.  Before  a  firing, 
the  surge  tank  and  motor  were  pre-pressurized  with  nitrogen.  The  tank 
provided  good  mean-pressure  regulation  during  the  test. 
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The  principal  measurement  in  the  system  was  the  pressure.  The  output  of 
the  transducer  system  (Kistler  Model  603A  transducer/Kistler  Model  504 
charge  amplifier)  was  recorded  simultaneously  by  a  magnetic  tape  recorder 
and  by  an  oscilloscope  camera.  A  secondary  measurement  was  the  thermo¬ 
couple  output.  This  signal  was  amplified  by  a  Kintel  DC  amplifier  and 
recorded  on  magnetic  tape  and  the  oscilloscope  as  well. 

The  experimental  system  was  operated  as  follows.  The  chamber  was  pressuri¬ 
zed  to  the  chosen  level.  Closing  an  igniter  switch  fired  the  igniters  and 
started  an  electronic  delay  generator  (thyratron  type).  This  delay  circuit 
fired  the  first  pulse-generating  squib  about  midway  through  the  firing 
by  closing  a  relay.  Near  the  end  of  the  test,  burnout  was  sensed  by  the 
thermocouple.  The  amplified  thermocouple  output  triggered  a  second 
electronic  delay  generator  which  fired  the  second  squib  about  90  milli¬ 
seconds  after  the  equilibrium  thermocouple  output  was  reached. 

1 

Three  kinds  of  squibs  were  used  in  this  program.  New  squibs  would  have 
been  prohibitively  expensive;  however,  several  kinds  of  squibs  were  avail¬ 
able  from  supply  at  no  cost.  The  first  two  kinds  were  used  until  the 
available  supply  was  exhausted.  Through  test  40,  a  Rocketdyne  initiator 
designated  as  NA5-26523-1A  was  used,  when  applicable.  For  tests  41  through 
140,  a  Holex  Corp.  squib  with  Holex  designation  of  2374  was  used,  when 
applicable.  For  tests  141  through  183,  a  second  Holex  squib  designated 
as  1785  was  used.  All  three  squibs  were  able  to  excite  oscillations  in 
the  motor.  The  first  squib  was  rather  unsatisfactory  because  it  intro¬ 
duced  excessive  damping,  apparently  due  to  condensed-phase  combustion  pro¬ 
ducts.  The  later  two  were  satisfactory  in  this  regard  and  seemed  to  be 
equivalent  for  purposes  of  this  experiment. 
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Several  CTPB/AP-type  nropellants  were  used  in  this  study.  They  are 
tabulated  in  Table  4.  Aluminum  content  was  varied  with  the  proportions 
of  other  constituents  remaining  constant.  The  aluminum  was  Valley 
Metallurgical  H-30  with  a  mass-mean  diameter  of  35  M  according  to  micro- 
merograph  analysis.  Composition  of  these  propellants  was  chosen  to  cover 
an  adequate  range  of  aluminum  loadings  and  to  be  similar  to  propellants 
which  have  been  studied  extensively  by  other  investigators. 


TABLE  4 

COMPOSITION  OF  PROPELLANTS  USED 


Composition  (weight  percent) 

Propellant 

A 

B 

c 

D 

F 

Binder 

20 

19 

17 

20 

19 

Fine  AP 

(20  fx,  nominal) 

39 

37 

33.2 

40 

38 

Coarse  AP 
(150  IX.  nominal) 

39 

37 

33.2 

40 

38 

Copper  Chromite 

2 

1.9 

1.7 

0 

0 

Aluminum 

0 

5 

15 

0 

5 
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EXPERIMENTAL  RESULTS 

The  experimental  pressure  data  were  interpreted  in  terms  of  a  growth 
constant,  O',  where  Ot  is  defined  by  the  time  dependence  of  the  data.  This 
dependence,  as  long  as  the  amplitude  is  not  too  large,  is  given  by  e0* 
cos  Wt.  At  high  amplitudes,  a  is  no  longer  constant  and  becomes  depend¬ 
ent,  on  the  amplitude.  If  the  oscillations  decay,  Of  is  negative,  (irowth 
constants  were  calculated  from  the  pressure  data  by  plotting  the  double- 
amplitude  (peak-to-peak)  against  time  on  semi logarithmic  coordinators. 

On  such  a  plot,  the  slope  of  the  dat3  represents  Of,  The  double-amplitude 
was  measured  from  oscillograph  traces  made,  at  reduced  speed,  from  the 
originally  recorded  data.  A  smooth  curve  was  drawn  through  the  oscilla¬ 
tory  peaks  and  this  was  taken  to  represent  the  amplitude  at  each  time 
increment. 

For  the  most  part,  the  data  were  well  represented  by  the  exponential  type 
of  time  dependence,  even  for  those  runs  in  which  the  oscillations  decayed 
very  rapidly.  In  some  of  the  low-frequency  tests,  the  waveform  contained 
significant  harmonic  distortion.  The  data  from  these  runs  were  closely 
low-pass  filtered  to  minimize  the  harmonic  content.  Scatter  in  the  com¬ 
pleted  values  of  O'  increased  with  frequency,  undoubtedly  much  of  this  can 
be  associated  with  the  very  high  decay  rates  encountered  in  this  region. 

The  observed  value  of  d  is  the  sum  of  several  contributions  (Ref.  26). 
Individual  contributions  can  be  evaluated  by  appropriate  subtraction  of 
the  a's  which  are  measured  in  different  kinds  of  measurements.  Such  an 
approach  is  used  to  evaluate  the  real  part  of  the  response  function, 

Re  (jp),  from  the  growth  constant  measurements.  The  relationship  is  (see, 
for  example,  Ref.  38,  39): 

»•  <?>  -  %  72  K  -  (®> 

s  2c 
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In  Eq.  T2the  assumptions  were  made  that  (l)  0!^  is  the  sum  of  a  contribution 
from  the  burning  surface  and  the  contributions  from  all  the  damping  pro¬ 
cesses,  and  (2)  tint  0^  represents  all  of  the  damping  contributions. 

The  additive  character  of  the  tt's  can  be  used  to  obtain  the  attenuation 
coefficient  for  the  increase  in  volumetric  damping  due  to  the  aluminum. 

This  damping  is  thought  to  be  essentially  entirely  cue  to  particulate 
damping.  The  quantity  O  ^  is  the  sum  of  a  contribution  due  to  the  pres¬ 
ence  of  aluminum  and  a  contribution  from  other  damping  processes.  The 
tt  from  these  "other  damping  processes"  was  assumed  equal  to  the  ob¬ 
tained  from  a  squib  test  with  no  aluminum  present.  Through  these  assump¬ 
tions  the  growth-constant,  data  were  evaluated  to  determine  the  changes 
in  response  function  and  volumetric  damping  for  the  aluminized  propellants 
relative  to  the  non-aluminized  propellants. 

Tests  were  made  at  mean-pressure  levels  of  200,  400,  and  600  psi  with 
motor  lengths  of  9,  1 5 ,  21,  and  41  inches  and  with  propellants  A,  D,  C, 

1),  and  I’  (Table  IV).  To  thoroughly  determine  the  growth  constants  for 
each  of  these  conditions,  a  prohibitively  large  number  of  tests  wou±d  be 
required,  ltecause  of  this,  emphasis  was  placed  on  propellants  A  and  C 
and  pressure  levels  of  200  and  600  psi  with  all  motor  lengths.  Fewer 
tests  were  made  at  the  other  conditions,  however,  the  results  of  these 
tests  are  entirely  commensurate  with  those  of  the  more  thoroughly  studied 
conditions.  In  addition  to  the  propellant  tests,  a  number  of  tests  were 
made  with  only  the  igniter  and  one  squib  being  fired.  These  tests  were 
made  to  evaluate  the  damping  introduced  by  the  squib,  the  igniter  being 
burned  to  heat  the  gas  in  the  motor  before  pulsing.  Tire  growth  constants 
obtained  from  all  tests  are  tabulated  in  Appendix  C.  Most  of  the  data 
are  plotted  in  Fig.  5  through  18. 
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On  Fig.  5  are  Pi  otted  the  growth  constants  obtained  from  ^  propellant 
at  200  psi.  In  this  case  the  motor  is  unstable  and  oscillates  throughout 
the  tests;  therefore,  is  a  positive  number.  The  curves  through  the 
data  were  used  to  evaluate  Re  (^)  for  this  propellant.  On  Fig.  6  are 
shown  similar  values  for  C  propellant  at  200  psi.  The  solid  curves  were 
drawn  to  represent  the  data.  In  this  case,  the  motor  is  stable  and  puls¬ 
ing  was  necessary.  Therefore,  the  plotted  values  include  a  damping  con¬ 
tribution  from  the  squib.  These  curves  may  be  used  to  calculate  lie  (^-) 
for  this  propellant.  Note  that  in  Eq.  22  the  coefficient  of  the  term 
(g^  -  flu '  )  is  approximately  the  same  for  each  propellant  since  the  burning 

rates  appear  equal  from  the  data.  If  the  lie  (^)  was  unaltered  by  the 
presence  of  aluminum,  as  found  by  Horton  (Ref.  6)  with  lightly  aluminized 
propellants,  the  vertical  distance  between  the  solid  curves  on  this  plot 
would  equal  the  difference  in  growth  constants  for  A  propellant  shown  on 
Fig.  5*  This  difference  from  Fig.  5  has  been  added  to  the  0^  curve  on 
Fig.  6  and  is  shown  as  the  dashed  curve.  The  dashed  curve  is  the  0!  ^ 
curve  that  would  be  obtained  if  lie  (^)  were  the  same  as  that  for  A  propellant. 

The  data  indicate  that  the  Re  (^)  Is  approximately  doubled  by  the  addition 
of  aluminum.  The  indication  that  the  aluminum  destabilizes  the  burning 
surfac  is  an  important  consideration  which  appears  in  all  of  the  experi¬ 
mental  results.  Refinement  of  the  experimental  technique  should  be 
carried  out  to  eliminate  all  possibility  of  systematic  errors  and  to 
further  verify  this  apparent  amplification. 


If  the  surface  has  been  greatly  stabilized  by  the  aluminum,  the  difference 
in  growth  constants  would  have  been  negative,  i.e.,  the  Cl  curve  would 
have  been  below  the  curve.  Clearly,  within  the  limitations  of  the 
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data,  the  surface  is  not  stabilized  by  the  aluminum.  Also  it  is  clear 
that  increased  overall  stability  of  aluminized  propellants  must  be  due 
to  increased  volumetric  damping.  It  is  likely  that  the  increased  volumetric- 
damping  is  principally  particulate  damping. 

The  increase  in  volumetric  damping  due  to  the  aluminum  was  calculated  as 
follows.  The  damping  associated  with  the  motor  and  squib  was  evaluated 
from  the  tests  with  only  a  squib  and  igniter  being  fired.  These  data 
ure  shown  in  Fig.  7-  At  the  two  higher  pressures  the  data  from  the 
A  propellant  tests  were  also  used.  All  of  these  data  are  described  by 
a  similar  curve  (Fig. H  and  14)so  that  these  additional  data  provide  some 
support  for  the  extrapolation  to  higher  frequencies  shown  in  Fig.  7. 

The  particulate  damping  was  calculated  from  the  Ot  data  by  subtracting 
the  extraneous  contribution  defined  by  the  curve  on  Fig.  7.  These  cal¬ 
culated  points  are  shown  on  Fig.  8.  It  is  apparent  that  the  particulate 
damping  was  substantial. 


Also  shown  on  Fig.  8  are  two  curves.  These  curves  were  calculated  from 
particulate-damping  theory  (either  the  Temkin  or  Epstein-Carhart  model) 
for  the  indicated  mean  particle  size.  The  important  observation  here  is 
that  a  fixed  mean-particle  size  will  not  describe  the  data  from  the  sev¬ 
eral  motors.  The  data  fall  in  the  region  where  the  mean  size  is  appro¬ 
priate.  The  data  indicate  that  the  size  distribution  changes  with  the 
length  of  the  motor,  llorton  (Ref.  7)  also  considered  a  change  in  particle 
size  to  be  likely. 

The  results  from  A  and  C  propellants  at  600  psi  were  much  the  same  as 
those  just  discussed.  On  Fig.  9  are  shown  the  growth  constants  obtained 
from  A  propellant.  In  this  case  the  motor  was  stable  and  pulsing  was 
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necessary.  Figure  10  shows  the  corresponding  results  for  C  propellant. 
On  this  plot  the  dashed  curve  tgain  represents  the  tt  ^  curve  that  would 
be  obtained  if  the  Re  (^)  were  the  same  as  for  A  propellant.  Clearly, 
the  results  indicate  the  aluminum  greatly  increases  Re  (jp) ,  even  more  so 
than  at  20G  psi.  The  damping  attributed  the  motor  and  squib  is  shown 
in  Fig.  11. 


The  results  of  the  particulate-damping  calculations  are  shown  on  Fig.  12 
On  this  plot,  the  circles  were  obtained  in  the  same  manner  as  the  values 
shown  on  Fig.  8.  The  second  set  of  points,  the  triangles,  were  calculated 
by  assuming  the  Re  (jp)  to  be  unchanged  by  the  aluminum,  i.e.,  in  the  same 
manner  as  that  used  by  Horton  (lief.  6).  Specifically,  the  decay  rate  de¬ 
fined  by  the  curve  for  A  propellant  was  subtracted  from  the  values 
for  C  propellant  to  give  the  particulate-damping  coefficient.  Also  shown 
on  this  plot  are  three  curves  calculated  from  particulate-damping  theory. 

As  before,  the  experimentally  determined  particulate  damping  is  not  de¬ 
scribed  by  a  single  particle-size  distribution.  Again  it  is  concluded 
that  the  size  distribution  changes  w.„h  motor  length. 

Results  for  A  propellant  at  400  psi  are  shown  on  Fig.  13.  The  damping 
attributed  to  the  squib  and  motor  is  shown  on  Fig.  14.  The  results  ob¬ 
tained  from  B  propellant  at  200,  400,  and  600  psi  are  shown  on  Fig.  15, 
l6,  and  17,  respectively.  All  of  these  results  are  very  similar  to  those 
discussed  previously.  Further  growth  constant  data,  which  have  not  been 
plotted,  are  tabulated  in  Appendix  C.  These  data  compare  well  with  those 
that  have  been  plotted. 


o 


Some  additional  observations  concerning  particulate  damping  may  be  made. 
Lxperimentally ,  the  effects  of  pressure  on  particulate  damping  appear 
small,  as  is  indicated  by  a  comparison  of  Fig.  8  and  12.  Particulate 
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FREQUENCY,  KILOCYCLES 

Comparison  of  Experimentally  Determined 
Particulate  Damping  With  Constant  Particle 
Size  Curves  for  C -Propellant  at  600  C 


Figure  14.  Decay-Rate  Data  for  l'-4  Squibs  at  400  psi 
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Figure  l1).  Decay-Rate  Data  for  B-Pro  pell  ant  at  'J(H)  psi 
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damping  theory  predicts  little  effect  of  pressure  on  damping  for  a  fixed 
particulate  size  (Ref.  34).  If  this  is  accepted  then  the  experimental 
results  cited  indicate  the  effect  of  pressure  on  particulate  size  to  he 
smal 1 . 

The  effects  of  particle  concentration  on  particulate  damping,  as  indicated 
by  the  experimental  results  obtained  are  approximately  in  accord  with 
theory.  On  1 ig.  18  is  shown  the  experimental  dependence  of  damping  on 
the  mass  ratio,  C  .  The  shape  of  the  curve  through  the  experimental 
val  ues  was  calculated  according  to  the  Xemlcin  model,  liq.  5.  At  other 
frequencies  and  pressures,  comparable  results  were  obtained. 

It  is  evident  that  the  pulsed  T-burner  employed  in  these  studies,  while 
subject  to  certain  limitations,  does  provide  a  powerf  1  tool  for  investiga¬ 
ting  the  effects  which  form  the  subject  of  this  report.  In  particular, 
the  pulse  concept  permits  the  desired  separation  of  surface  and  chamber 
volume  effects,  allowing  this  to  be  accomplished  at  aluminum  loadings  in 
the  range  employed  in  conventional  propellants  and  inaccessible  to  con¬ 
ventional  T-burner  investigation.  The  results  of  this  separation  show 
that  the  dominant  mechanism  by  which  aluminum  operates  to  promote  stability 
is  particulate  damping  in  the  gas  phase.  They  further  suggest  that  the 
aluminum  may  actually  tend  to  destabilize  the  burning  surface.  Although 
for  the  motors  considered  in  this  report  the  damping  influence  outweighs 
the  destabilizing  tendency,  conceivably  in  motors  possessing  a  higher 
ratio  of  burning  surface  area  to  chamber  volume  the  overall  stabilizing 
effect  of  aluminum  could  be  considerably  reduced. 
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CONCLUSIONS  AND  RDC OMMFNDAT IONS 

The  investigation  reported  has  been  considerable  success  in  throwing 
light  upon  the  role  of  aluminum  in  combustion  instability  of  solid  pro¬ 
pellants.  This  effort  represents  the  first  successful  attempt  to  obtain 
quantitative  information  on  the  stability  characteristics  of  propellants 
which  contain  sufficient  aluminum  to  be  of  practical  interest.  Farther, 
the  broad  scope  of  the  program  has  been  fruitful  in  developing  interest¬ 
ing  areas  for  further  research. 

The  analytical  phase  has  indicated  that: 

1.  The  dominant  stabilizing  influence  associated  with  aluminum  is 
due  to  particulate  damping  in  the  chamber  volume. 

2.  Various  kinds  of  oscillatory  motion  of  the  droplet  are  not  likely 
to  play  an  important  role. 

3.  Drop  shattering  is  probably  unimportant  unless  the  shattering 
frequency  is  greater  than  one  would  reasonably  estimate.  Photo¬ 
graphic  studies  are  needed  to  determine  the  extent  of  shattering. 

A.  The  effect  0?  aluminum  in  modifying  the  interaction  between  the 
acoustic  field  and  the  combustion  may  be  significant  (though 
probably  still  small  with  respect  to  particulate  damping).  How¬ 
ever,  analysis  of  this  interaction  ras  found  to  be  completely 
intractable,  since  it  would  require  analysis  of  the  complete 
heterogeneous  combustion  problem. 

5.  Nonacoustic  instability  is  closely  related  to  acoustic  instability. 
This  relationship  suggests  the  need  to  apply  the  current  knowledge 
of  the  latter  to  improve  understanding  of  the  former. 
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f>.  The  Den  is  on -Baum  cora.mfction  model  can  be  used  to  quantitatively 
describe  both  ac  mstic  and  nonacoustic  instability  data.  The 
nonacoustic  instability  data  included  data  from  10  percent 
aluminized  propellants.  Farther  use  of  this  model  to  explain 
experimental  observations  is  suggested. 

7.  i^uasi-periodic  shedding  of  droplets  from  the  burning  surface  is 
likely  to  occur  with  a  frequency  near  the  maximum  in  the  response 
function  curve.  Experimental  consideration  of  this  relationship 
is  needed. 

During  the  experimental  study,  instability  data  were  obtained  from  pro¬ 
pellants  containing  up  to  15  weight  percent  aluminum  and  at  pressure  levels 
as  high  as  600  psi.  Although  the  data  arc  of  limited  precision,  important 
information  was  obtained.  The  data  indicate  that: 

1.  Particulate  damping  is  the  domonant  mechanism  whereby  aluminum 
stabilizes  acoustic  instability  at  all  aluminum  loadings.  Sub¬ 
stantial  particulate  damping  was  found, being  measurable  at 
frequencies  as  low  as  200  cps, 

2.  The  inclusion  of  aluminum  destabilizes  the  burning  surface  over 
the  entire  frequency  range;  i.e.,  the  surface  has  increased 
ability  to  amplify  a  pressure  disturbance. 

5.  The  pa»-ticle-Bize  distribution  varies  with  T-burner  length. 

4.  The  pressure  effect  on  particulate  damping  is  small,  which  im¬ 
plies  that  the  pressure  dependence  of  the  particle- -»e  dis¬ 
tribution  is  small. 

5.  The  effect  of  aluminum  concentration  is  reasonably  in  accord 
with  theory. 


54 


R-6654 


«oc«E'rE>\'roie 


A  DIVISION  OF  NORTH  AMERICAN  AVIATION 


I  N  ( 


s* 


Further,  more  precise  experimental  work  is  needed  to  fully  verify  each 
of  these  indications.  The  destabilization  of  the  burning  surface  is 
particularly  important,  both  for  acoustic  and  nonacoustic  instability. 

Finally,  the  experimental  study  has  demonstrated  the  effectiveness  of  a 
pulsed  T-burner  as  a  means  of  studying  propellants  too  stable  to  drive 
the  T-burner. 
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APPENDIX  A 


NOMENCLATURE 
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m 

C 

P 
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l 

■*  X 


D.  . 
ij 


part 


"them 


defined  by  Eq.  D-ll 

isentropic  sound  velocity,  yp Q/ PQ 

ratio  of  total  particulate  mass  per  unit  volume  to  gas  density 

heat  capacity  at  constant  pressure  for  the  gas 

heat  capacity  of  particulate  matter 

heat  capacity  (average)  of  solid  propellant 

isentropic  sound  velocity,  yp Q/PQ 

isentropic  sound  velocity  in  droplet 

characteristic  exit  velocity  of  the  gas  (isentropic) 

motor  port  diameter 

mean  particle  size,  defined  by  Eq.  10 
diameter  of  droplet  located  on  burning  surface 
diameter  of  droplet  located  in  main  gas  volume 
activation  energy  for  overall  gas-phase  reaction 
activation  energy  for  surface  vaporization 
defined  by  Eq.  D-ll 
rate  of  energy  dissipation 
E  due  to  particulate  damping 
thermal  contribution  to  E 
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n 

Pr 

Po 

P 

A 
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viscous  contribution  to  E 
oscillatory  frequency 
imaginary  part  of  enclosed  quantity 
(-1)1/2 

ratio  of  burning  surface  area  throat  area 
thermal  conductivity 

k  of  droplet 
motor  length 

ratio  of  chamber  volume  to  throat  area 

motor  length 

time-average  mass  flux 

perturbed  mass  flux 

number  of  drops  in  burning  surface 

number  of  drops  in  chamber  volume 

burning  rate  exponent 

integer 

distribution  in  particle  size 
M  Cp/k 

time  average  pressure 
perturbed  pressure 
amplitude  of  p 
defined  by  Eq.  B-ll 
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r>. 


a 


nb 


a 


ideal  gas  constant 

real  part  of  enclosed  quantity 

burning  rate  of  propellant 

radius  of  particle 

mean  particle  radius 

defined  by  Eq.  B-ll 

burning  surface  area 

time  average  temperature 

perturbed  temperature 

Tq  cf  combustion  products 

initial  temperature  of  unburned  propellant 

Tq  at  burning  surface 

T  of  combustion  products 

time 

time-average  gas  velocity 
chamber  volume  (gas  phase) 
growth  constant 
thermal  diffusivi ty 

growth  constant  obtained  durng  burning 
thermal  expansion  coefficient  of  droplet 
growth  constant  obtained  after  burn-out 
reduced  attenuation  coef f  icient,  Eq.  1 
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=  (c-  j  a 

=  compressibility  of  droplet 

-  ratio  of  heat  capacities  for  gas  phase,  C  /Cy 

-  ratio  of  heat  capacities  for  droplet 
=  propagation  constant,  Eq.  B-5 

=  ratio  of  gas  density  to  particle  density 


phase  angle 
defined  by  Eq,  B-ll 


=  gas  viscosity 


droplet  viscosity 


-  3.14159 
=  density 

time-average  p 
=  time-average  gas  density 
tine-average  droplet  density 
=  density  (average)  of  solid  propellant 


surface  tension 


=  defined  by  Eq.  11-11 

(d;  pj)  (8m) 
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(-,  ft-  r  y  (2  cp) 

virtual  acoustic  admittance  of  nozzle 
^  viscous  dissipation  function 
=  angular  frequency,  2M 
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APPENDIX  B 

ANALYSIS  OF  NONACOUSTIC  INSTABILITY 

Quas i-periodic  agglomeration  and  shedding  of  aluminum  from  the  burning 
surface  at  lov  frequencies  has  often  been  reported  (Ref.  20  and  21 ).  Also 
it  has  often  been  suggested  that  this  process  may  act  to  drive  nonacoustic 
instability.  In  the  nonacoustic  frequency  range  particulate  damping  is 
expected  to  be  unimportant  due  to  its  strong  dependence  on  frequency. 
Therefore,  the  aluminum  may  have  a  destabilizing  influence  in  this  fre¬ 
quency  region.  The  effect  of  the  aluminum  is  difficult  to  estimate 
analytically.  Some  progress  was  made  in  determining  the  source  of  the 
quasi-periodicity,  however.  In  doing  so,  a  detailed  analysis  of  non¬ 
acoustic  instability  was  developed.  Significantly,  it  was  found  that 
nonacoustic  instability  is  intimately  related  to  acoustic  instability 
rather  than  completely  divorced  as  other  investigators  have  concluded. 

The  dynamic  behavior  of  a  motor  exhibiting  nonacoustic  instability  was 
analyzed.  Previous  analyses  of  nonacoustic  instability  have  made  use 
of  the  ballistic  equation,  in  perturbation  form,  and  a  model  for  the 
combustion  behavior  (Ref.  20  and  21) .  Following  this  approach,  but  without 
specifying  the  combustion  model,  a  characteristic  equation  was  obtained 
for  the  frequency. 

It  is  assumed  that:  (l)  the  mode  is  a  bulk  mode;  (2)  the  gas  behavior  is 
isentropic;  the  isothermal  case  is  obtained  by  setting  y  to  unity;  (3)  the 
isentropic  expression  for  flow  through  the  nozzle  is  valid  in  the  per¬ 
turbed  case.  This  is  equivalent  to  the  zero  length-nozzle  approximation 
for  the  acoustic  admittance  of  a  nozzle.  The  approximation  is  good  be¬ 
cause  the  equivalent  wavelength  is  infinite;  and  (4)  the  rate  of  volume 
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increase,  due  to  time-average  burning  is  negligible.  With  these  assump¬ 
tions  the  perturbed  equation  becomes: 


L*  c*  , 
o  d 


2  dt  V  p 
c  r0 


(B-l) 


Assuming  an  exponential  time  dependence,  i.e.,  eJP  .  Equation  B-l  is 
rearranged  to  give 


~  -  L*  c*  *  i  x  ~ 

—  -  =  (  j  fj - -JL.  +  Jtl  ^  £- 

m  v  ^  2  2  y  y  v 


(B-2) 


If  use  is  made  of  the  definition  of  the  response  function,  in  perturbation 
form,  we  obtain: 


L*  c* 


(B-3) 


Therefore,  the  response  function  explicitly  defines  the  time  dependence, 
R.  The  frequency  of  oscillation  is  determined  by  the  imaginary  part. 

Im  {tt'c},  and  whether  the  oscillations  grow  or  decay  is  determined  by 
the  real  part,  Re  {j|/c}. 

Equation  B-3  fan  also  be  obtained  by  an  acoustic  approach.  For  low- 
velocity  mean  flow  in  the  main  cavity,  the  mean  flow  may  be  neglected 
as  long  as  it  is  taken  into  account  at  the  boundaries  by  using  virtual 
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acoustic-admittance  values  at  these  positions.  Assuming  the  wave  behavior 
to  be  described  by  the  conventional  wave  equation,  solxitions  were  obtained 
for  two  motor  geometries. 


The  first  geometry  considered  was  cylindrical  with  propellant  being  burned 
on  the  curved  boundary,  a  rigid  boundary  at  the  head  end,  and  a  sonic 
nozzle  at  the  aft  end.  Without  difficulty  it  is  found  that 


Y —  4. — u  4  a,  i  s  _  v 
p  $  J'  '  *n 


(3-4) 


where  y'  is  given  by 


,A.  (£.v2 


£  -  y2 
2  r 
.c  / 


j 1/2  I 


<£  ,2 


The  argument  of  the  Bessel  functions  is  normally  sufficiently  small  that 
Eq.  B-5  becomes: 


<  £sl  a,  (£ .  v- 


o  r  4/9  \  2 
s  '  c 


(B-6) 


This  equation  may  be  solved  for  y'  ,  i.e., 


i  &  ^  •  c  a 

V  c  J  2L*  c  *  € 
0 


(B-7) 
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Equation  B-4  conjunction  with  Eq.  B-7,  determines  the 

complex  frequency,  P. 

In  the  case  of  interest,  the  virtual  nozzle  admittance  is  sufficiently 
small  that 


tanh  j  y'  ■£•  ~  J  (y’^  ~  nil)  n  =  0,  1,2,  ... 


(B-8) 


The  values  of  n-  correspond  to  various  longitudinal  modes.  With  the  zero- 
length  nozzle  equation,  and  setting  n  =  0  for  the  zeroth  mode,  it  is 
found  from  Eq.  3-8  that 


V  L*  c  * 


c 


(B~9) 


which  is  identical  with  Eq.B-3  .  Thus  the  zeroth  mode  corresponds  to 
the  bulk  mode  previously  treated  with  the  ballistic  equation.  Note  that 
for  a  closed  cavity  with  no  gains  or  losses,  /?  =  0  for  the  zeroth  mode. 
Therefore,  there  is  no  analogous  mode  in  the  simple  case.  Also,  it  is 
observed  that  several  approximations  were  made  in  order  to  obtain  the 
results  given  in  Eq.  B-9.  All  of  the  approximations  are  usually  valid; 
however,  if  any  of  the  approximations  are  invalid  then  the  more  general 
equations  must  be  used.  Indeed,  in  the  latter  case,  use  of  the  perturbed 
ballistic  equation  is  incorrect. 
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The  3eeond  configuration  considered  vas  also  a  cylindrical  geometry  but 
the  curved  walls  were  rigid,  propellant  was  burned  at  the  head  end  and 
a  sonic  nozzle  was  located  at  the  aft  end.  Again  using  the  zero-length 
nozzle  approximation,  solution  of  the  acoustic  problem  leads  to  the 
characteristic  equation 


(B-10) 


With  reasonable  values  of  Kq>  the  zeroth  mode  is  again  given  by  an 
equation  identical  with  Eq.  B- 3. 


The  foregoing  results  emphasize  the  relationship  between  acoustic  and 
nonacoustic  instability.  Indeed,  they  suggest  that  combustion  models 
for  acoustic  instability  (e.g.Ref.  24  and  25)  may  apply  to  nonacoustic 
instability  as  well.  The  value  of  this  suggestion  was  tpcted  with  the 
Denison-Baum  model  for  acoustic  instability,  (Ref.  251 • 


For  an  exponential  time  dependence  (i.e.,  e^”  ),  the  Denison-Baum  model 
may  be  solved  to  obtain  the  following  response  function: 


-n  A  A _ 

q  +  s  "  (l  +  s) 


(B-ll) 


where 


S  =  1/2  [l  +  (1  +  4 j Ar) 1/2 j 


q  =  1  +  A  (l  -  A) 


A  ► 
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'  £  C  (T  -  T.J 

-  S  WO  1 


£  =  n  +  1  +  Ef/2HTfo 


E  y  T.  N 

A  «  _2L.fl  _  -i_> 
RT  V1  T  J 
wo  wo 


T-  =  flame  temperature 


T.  =  bulk  temperature  of 
propellant 


The  parameters  X,  n,  and  S  are  equivalent  to  Oc,  n/2,  and  X,  respectively, 
in  Denison-Baum  nomenclature. 

Hart  and  McClure  (Ref,  2 6)  indicate  that  the  Denison-Baum  results  are 
quite  similar  to  their  own.  However,  no  direct  comparison  of  either 
model  has  been  made  with  experimental  data.  Therefore,  extensive  calcu¬ 
lations  were  made  with  Eq.  B-ll  to  test  its  ability  to  describe  response- 
function  data. 

Perhaps  the  best  data  available  from  which  to  calculate  the  Re  j— are 
reported  by  Horton  in  Ref.  27  from  T-bumer  firings.  Using  Horton's 
data,  the  Re  was  calculated  from  the  equation 


Ee  f-4  = 


2r  p  c 

's 


a,  -  a 

b  n. 


(B~12) 


Similar  equations  are  widely  used;  an  equivalent  equation  is  developed 
by  Cantrell  (Ref.  18).  The  data  for  were  plotted  as  a  function  of 
frequency  and  a  smooth  curve  drawn  through  the  data.  An  example  cf  this 
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procedure  is  shown  in  Fig.  B-l*  In  this  figure,  lines  are  drawn 
connecting  pairs  of  data  points;  each  connected  pair  of  points  represents 
the  initial  and  final  values  of  0^  reported  by  Horton.  The  smooth  curve 
through  the  data  has  been  drawn  to  approach  zero  as  the  frequency  approaches 
zero.  This  was  done  because,  whatever  the  damping  mechanism,  the  damping 
coefficient  is  expected  to  vanish  at  zero  frequency.  Th<  smooth  curve  was 
then  used  to  define  a  value  of  to  be  used  in  Eq.B-12  with  the  Qb  data, 
thereby  correcting  for  frequency  shifts  between  the  %  and  measurements 
as  well  as  smoothing  the  data. 

Thus,  experimental  values  of  Re  were  calculated  for  four  different 

propellants  and  are  presented  in  Fig.  B-2  through  B-5  •  The  data  were 
obtained  at  a  mean  pressure  of  200  psi.  The  propellants  were  designated 
as  A-13,  A-14 ,  A-15  and  A-16  by  NOTS.  Each  was  of  the  PBAA/AP  type. 

A-15  and  A-16  contained  copper  chromite  catalyst  while  A-13  and  A-14 
contained  none.  Propellants  A-13  and  A-15  contained  76-pcrcent,  80-micron 
AP,  while  A-14  and  A-16  contained  the  same  amount  of  15-micron  AP.  On 
these  plots  are  shown  curves  through  the  experimental  points  which  were 
calculated  from  the  theoretical  expression.  The  curves,  which  describe 
the  experimental  values  rather  well,  were  obtained  by  varying  the  parame¬ 
ters  A  <  d  q  in  Eq.  B-il  until  the  best  fit  was  obtained.  Although  a 
rather  good  description  of  the  experimental  values  is  obtained,  the  model 
is  somewhat  unsatisf actory  since  the  curves  correspond  to  unusually  high 
activation  energies  for  the  surface  decomposition  (65  to  100  kcal/mole-K) 
However,  this  activation  energy  must  be  considered  to  be  an  effective  value 
„  _  since  the  model  approximates  the  heterogeneous  propellant  by  a  homogeneous 

solid.  Therefore,  perhaps  the  activation  energy  required  to  fit  the 
experimental  Re  values  is  not  unreasonable .  Also,  there  is  no  doubt 

that  the  model  can  qualitatively  describe  the  experimental  values.  It 
was  concluded  that  the  model  was  very  successful  for  these  purposes. 

*  v 
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It  should  be  mentioned  that  the  high-frequency  limitation  of  the  Denison- 
Baum  model  is  no  problem  for  the  propellants  considered.  According  to 
the  criterion  specified  by  Hart  (Ref.  26),  the  model  was  estimated  to 
become  invalid  above  Ur  ~  240  for  each  propellant.  Nearly  all  of  the 
data  fall  below  this  value  of  U5T. 

A  second  detailed  comparison  was  made  with  nonacoustic  instability  data 
of  Beckstead  (Ref.  22).  He  has  studied  the  behavior  of  several  propellants 
in  a  low  L*  motor.  In  the  motcr,  propellant  was  burned  in  one  end  of  a 

cylindrical  cavity  c-nd  a  sonic  nozzle  was  contained  in  the  other  end.  He 

reports  values  of  L*,  frequency,  and  mean  pressure  during  instability, 
and.  therefore,  provides  sufficient  information  to  calculate  Im  J  — j- 
from  Eq.  i>-3  .  Other  investigators  have  not  reported  these  data  in 
sufficient  detail  to  caltulate  Im  In  Fig.  R-6  are  presented  the 

Im  (~j  obtained  from  Beckstead 's  results  for  three  of  his  propellants 
and  plotted  in  the  format  used  by  him.  In  this  format  a  mean-pressure 
effect  in  the  data  is  apparently  suppressed.  Beckstead 's  propellants 
wer*  of  the  PBAA/AP  type  with  two  percent  copper  chromite  catalyst. 
Propellants  F,  TF,  and  XF  contained  0.  5,  and  10  percent  aluminum  and 

80,  75.  and  70  percent  AP,  respectively.  On  Fig.  B-6  is  shown  a  curve 

through  the  experimental  values  calculated  from  Eq.  3^3  .  Again  the 
values  of  A  and  q  in  Eq.  B-ll  were  varied  until  a  best  fit  was  obtained. 

As  with  the  T-burner  results,  the  model  describes  the  experimental  points 
rather  well.  However,  contrary  to  the  comparison  with  T-burner  results, 
an  activation  energy  for  the  surface  decomposition  as  low  as  27  kcal/mole-°K 
could  be  used  to  calculate  the  best-fit  curve.  Again  it  was  concluded 
that  the  model  describes  the  experimental  results  rather  well. 

Since  the  Denison-Baum  model  seemed  fully  adequate  to  at  least  qualitatively 
describe  unstable  combustion,  several  qualitative  features  observed  in  non¬ 
acoustic  instability  studies  were  considered.  In  each  case  an  adequate 
explanation  was  obtained.  An  example  is  the  preferred  frequency  reported 
by  N0TS. 
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NOTS  uses  two  kinds  of  experimental  motors  in  their  preferred  frequency 
work.  One  is  the  L*  motor,  which  is  similar  to  Beckstead's  motor,  and 
the  other  a  very  long  T-burner,  denoted  as  the  LE  burner  (Ref.  20,  28, 
and  40.  In  the  L*  burner,  it  is  found  that  frequency  correlates  well 
with  mean  pressure  in  the  motor,  only  a  limited  band  of  frequencies 
beirg  observed.  Beckstead  has  pointed  out  that  there  is  an  L*  dependence 
in  this  band  of  frequencies  (Ref.  22  )•  Yount  (Ref.  41  )  has  reported  nc 
L*  dependence;  however,  this  may  be  attributed  tj  the  range  of  L*  values 
in  which  he  worked.  Yount's  data  also  correlate  on  a  frequency-mean 
pressure  plot.  Beckstead  found  a  similar  relationship. 

The  L*-(nonacoustic)  instability  data  may  be  explained  as  follows.  When 
the  L*  and/or  pressure  are  reduced,  the  frequency,  or  U)  r,  is  increased 
in  accord  with  Im  {u/e}  for  the  particular  propellant.  Simultaneously, 
the  Re  [ti/e]  is  increased  as  u)  T  increases.  When  Re  {14/e]  exceeds  the 
losses  in  the  system,  i.e.,  y+\/2y  (or  unity  for  an  isothermal  system 
if  the  only  loss  is  the  nozzle),  the  motor  becomes  unstable.  Continued 
increase  of  W  r  results  in  large  values  of  Rc  {p/?},  and  the  motor 
becomes  so  unstable  that  chuffing  results.  This  description  is  consistent 
with  experimental  observations. 

In  the  second  kind  of  experiment  carried  out  by  NOTS,  propellant  is 
burned  in  a  T-burner  ranging  in  length  up  to  60  feet.  In  this  system, 
the  mean  pressure  is  allowed  to  increase  slowly  during  the  run.  The 
frequency  is  determined  by  the  acoustics  of  the  cavity.  It  is  found 
that  instability  develops  when  the  pressure  becomes  sufficient  and  then 
ceases  at  a  higher  pressure  (Ref.  20  ,  28  ,  and  42). 

In  this  motor,  the  losses  are  determined  by  the  acoustics  of  the  cavity, 
these  being  frequency  depended .  Viewed  in  terms  of  the  Re  {i//e},  the 
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experiment  begins  at  a  large  value  of  l 0  7\  because  of  the  lot,-  pressure 
and  fixed  frequency.  As  the  pressure  rises,  lt)r  diminishes  and  Ho  lye} 
increases  tovard  the  maximum.  There  will  be  a  value  of  Re  {4/e}  which 
corresponds  to  the  losses  in  the  system;  when  this  value  is  exceeded,  the 
motor  will  oscillate.  As  U)T  continues  to  decrease,  the  oscillatory 
amplitude  will  increase  and  then  diminish.  When  Re  {u/ej  again  drops 
below  that  necessary  to  overcome  the  losses,  the  oscillation  will  cease. 

If  the  motor  length  is  changed,  the  frequency  is  changed,  and  also  the 
losses.  The  motor  will  again  oscillate  in  the  range  of  (Or  values  for 
which  He  lu/(}  exceeds  that  corresponding  to  the  losses.  If  the  losses 
are  sufficiently  large,  there  will  be  no  value  of  Re  {(j/f.)  at  which  the 
losses  are  exceeded. 

The  preferred  frequencies  will  coincide,  for  the  two  kinds  of  experiments 
only  if  the  losses  ure  equivalent  in  the  two  systems.  The  preferred  fre¬ 
quency  plot,  i.e.,  the  frequency-mean  pressure  plot,  shows  the  variation 
of  the  Wr  at  the  stability  limit  with  mean  pressure.  If  iu/e]  and  the 
losses  vary  only  weakly  with  mean  pressure,  the  preferred  frequency  will 
approximately  correspond  to  a  constant  WT  at  the  stability  limit. 

NOTS  has  also  reported  observing  a  phase  difference  between  the  oscillatory 
light  emission  from  the  propellant  surface  and  the  oscillatory  pressure 
(lief .  20  ,  28  .  29  ,  and  42).  This  phase  shift  has  been  explained  in  terms 
cr  the  preferred  frequency  concept.  The  observation  is  also  explainable 
in  terms  of  the  current  analytical  model. 

In  the  long  T-burner  NOTS  observed;  en  oscillating  light  signal  which 
lags  the  pressure  as  the  oscillations  begin;  as  the  pressure  oscillations 
reached  a  maximum,  the  phase  angle  became  zero;  and,  finally,  a  phase 
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load  was  observed  while  the  oscillation  diminished  (Ref.  28,  29,  and  42). 
In  the  L*  burner,  only  a  phase  lead  is  observed.  An  explanation  may  be 
obtained  if  it  is  assumed  that  the  light  oscillations  correspond  to 
oscillations  in  fleme  temperature.  The  Denison-Baum  model  gives: 


(<-“) 


fo 


45 
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JEL 

P„ 


(B-13) 


From  Eq.F—13  ,  the  phase  angle  between  the  perturbed  temperature  and 
pressure  is  given  by: 


tan  6  - 


\ 

Im 


(3-14) 


As  noted  previously,  in  the  long  T -burner,  OIT  is  large  initially  and 
diminishes.  In  relation  to  Eq.  B-14,  the  numerator  of  the  right  side 
is  negative  and  the  denominator  positive  at  large  Wr.  As  01  T  diminishes, 

Im  diminishes,  becoming  zero  at  the  01 T  corresponding  to  the  maxi¬ 

mum  in  Re  {u/c)  while  the  denominator  remains  positive.  Both  the  numera¬ 
tor  and  denominator  are  positive  for  values  of  01  T  less  than  that  of  the 
maximum  in  Re  {n/ej.  Thus,  Eq.  B-14  predicts  phase  angle  behavior  identical 
with  that  observed,  i.e.  ,  the  phase  angle  is  initially  negative,  then 
diminishes  and  becomes  positive. 

With  regard  to  the  low  L*  data,  the  motor  is  unstable  only  in  the  vicinity 
of  the  stability  limit,  i.e.,  in  the  region  where  Re  (u/cj  =  y+l/2y, 
assuming  nozzle  losses  dominate.  In  this  region,  the  phase  angle  is 
positive,  conesponding  to  the  observed  values. 
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At  N0T8  and  elsewhere,  quosi-periodic  shedding  of  aluminum  from  the 
burning  surface  of  aluminized  propellants  has  been  observed.  The  nominal 
frequency  of  this  shedding  is  in  the-  range  of  nonacoustic  instability. 

Such  shedding  is  observed  when  there  is  no  acoustic  cavity  in  the  vicinity 
of  the  propellant;  i.e.  ,  the  propellant  is  burned  in  the  open  or  in  a 
bomb.  However,  the  dense  cloud  of  burning  aluminum  leaving  the  burning 
surface  will  be  able  to  reflect  local  disturbances  generated  by  the  com¬ 
bustion  bach  to  the  surface  where  the  disturbance  can  be  amplified.  The 
frequency  components  of  the  disturbance  near  the  maximum  in  Re  [n/e] 
(actually  the  growth  constant  for  the  surface)  will  be  amplified  to  the 
greatest  degree.  It  seems  likely,  therefore,  that  the  frequency  of  the 
shedding  will  tend  toward  that  corresponding  to  the  maximum  in  the 
response  function. 

From  this  study  of  nonacoustic  instability,  it  was  concluded  that  a 
rather  good  analytical  description  of  the  phenomena  had  been  developed. 
Also,  the  utility  of  the  Denison-Baum  combustion  model  is  very  evident. 
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TABLE  C-I 

GROWTH  CONSTANT-FREQUENCY  DATA  OBTAINED  BY  FIRING  IGNITER 
ONLY  AND  PULSING  WITH  A  NO.  2374  SQUIB 


L, 

inches 


Test 

Nc. 

V 

psig 

‘  °nb’ 

seconds 

127 

200 

36.0 

128 

200 

46.0 

129 

200 

56.0 

96 

200 

19.0 

98 

200 

17.0 

92 

200 

13.0 

93 

200 

15.0 

94 

200 

13.0 

108 

200 

4.0 

109 

200 

3.3 

110 

200 

4.9 

130 

400 

39.0 

131 

400 

42.0 

132 

400 

47.0 

77 

400 

16.0 

79 

400 

13.0 

80 

400 

11.0 

44 

400 

16.0 

83 

400 

7.9 

84 

400 

8.2 

85 

400 

9.8 

111 

400 

5.6 

114 

400 

6.5 

115 

400 

6.2 

133 

600 

24.0 

134 

600 

28.0 
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TABLE  c_4 


GROWTH  CONSTANT-FREQUENCY  DATA  OBTAINED  FROM  C-PSOPELLANT 


‘V_1  fb’  ‘  “ '“I,  W 

seconds  cps  seconds  cps 


L, 

inches 


200 

230 

200 

310 

200 

350 

200 

350 

200 

180 

200 

140 

200 

200 

200 

150 

200 

180 

200 

30 

200 

36 

200 

41 

400 

110 

400 

140 

400 

a  yJyJ 

400 

82 

600 

240 

600 

210 

600 

260 

600 

260 

600 

180 

600 

140 

600 

160 

600 

120 

600 

56 

600 

61 

600 

42 

600 

40 
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TABLE  C-5 


GROWTH  C  ONSTAOT-FHEQUINCY  OBTAINED  FROM  D-PR0PFL1ANT 


Test 

No. 

V 

psig 

V 

seconds  * 

V 

cps 

-“nB 
seconds  ^ 

fnb’ 

cps 

L, 

inches 

75 

200 

10.0 

1150 

13.0 

R5 

15 

89 

200 

8.1 

1130 

16.0 

■ 

15 

200 

12.0 

li:0 

18.0 

980 

15 

74 

4.9 

810 

14.0 

21 

91 

3-1 

780 

9.4 

1 

21 

5.1 

310 

7.3 

41 

!  n, 

2.2 

280 

8.0 

111 

41 

TABLE  C-6 


GROWTH  C ONSTANT-FREQUQJCY  DATA  OBTAINED  FROM  E-PROPELLANT 


Test 

No. 

V 

psig 

-*b 

seconds 

fb’ 

cps 

-  a 

nb 

,  -1 
sec  onus 

fnb 

cps 

L, 

inches 

103 

200 

J  i 

1160 

67 

830 

15 

105 

200 

47 

1160 

67 

15 

112 

200 

28 

I18G 

91 

15 

102 

21 

820 

52 

750 

21 

104 

200 

32 

800 

61 

21 

106 

28 

810 

57 

740 

21 
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■  >  abstract  An  analytical  and  experimental  program  to  investigate  the  effects  of 
aluminum  on  solid-propellant  combustion  instability  -was  conducted.  The  ana¬ 
lytical  phase  vas  devoted  to  a  dplevmi nation  of  the  energy  dissipation  rate  ; 
associated  with  possible  damping  mechanisms.  The  relative  magnitude  of  these  \ 
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of  these  experiments,  particularly  with  respect  to  the  damping  mechanisms  con¬ 
sidered  in  the  analysis.  Some  information  applicable  to  nonacoustic  instability 
was  obtained  as  veil.  (U) 
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